The Luttinger model following a sudden interaction switch-on 
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The evolution of correlations in the exactly solvable Luttinger model (a model of interacting 
fermions in one dimension) after a sudden interaction switch-on is analytically studied. When the 
model is defined on a finite-size ring, zero-temperature correlations are periodic in time. However, 
in the thermodynamic limit, the system relaxes algebraically towards a stationary state which is 
well described, at least for some simple correlation functions, by the generalized Gibbs ensemble 
recently introduced by Rigol et al. cond-mat/0604476 . The critical exponent that characterizes 
the decay of the one-particle correlation function is difi^erent from the known equilibrium exponents. 
Experiments for which these results can be relevant are also discussed. 

PACS numbers: 



Experiments with cold atomic gases are motivating re- 
search into problems that, previously, would have looked 
highly academic. One such problem concerns the evolu- 
tion of a quantum many-body system where interactions 
(or other parameters of the system) are time-dependent. 
An example is an interaction quench: an experiment 
where the strength of interactions is suddenly changed. 
This type of experiment is nowadays feasible thanks to 
the phenomenon known as Feshbach resonance P, 
which allows to tune the strength and sign of interac- 
tions in a cold atomic gas by means of a magnetic field. 
If the applied magnetic field is time dependent, the inter- 
actions become time-dependent. Alternatively, in optical 
lattices , it is possible to change the lattice parameters 
in a time-dependent fashion, which effectively amounts to 
varying the ratio of the interaction to the kinetic energy 
in time. On the theory side, the recent development of 
extensions of the density-matrix renormalization group 
(DMRG) algorithm has spurred the interest 

in understanding the properties of quantum many-body 
systems out of equilibrium and, in particular, in the dy- 
namics following a quench. 

Because of these new possibilities, the evolution of ob- 
servables and correlations following a sudden change of 
the system parameters is attracting rnuch theoretical in- 
terest 11 EE [nl [11 [11 Hill [11 111 El One interest- 
ing question that has been raised by a recent experiment 
in an array of ID cold atomic gases ^1| is whether after 
a quench a system possessing an infinite number of inte- 
grals of motion can exhibit relaxation towards a steady 
state or not. This question has been analyzed by the au- 
thors of Ref. Ell J who have numerically shown that the 
steady state of an integrable gas of hard-core bosons is 
described by a generalized Gibbs distribution that max- 
imizes the entropy with all possible constraints imposed 
by the existence of the (infinite number of) integrals of 
motion. Here the effect of suddenly turning on the in- 
teractions in the Luttinger model is analytically studied. 
It is shown that, when the model is defined on a finite- 
size ring, the asymptotic form of the two-point one-body 
and density correlations at zero temperature is periodic 



in time, and therefore the system exhibits no relaxation 
to a steady state with time-independent properties. In 
the termodynamic limit, however, the same correlation 
functions relax to a steady state, whose properties are 
different from those of the ground state. Indeed, the de- 
cay of the one-particle correlations with distance is gov- 
erned by a critical exponent which is different from the 
known equilibrium exponents. Interestingly, one-particle 
and density correlations in the steady state can be ob- 
tained using the generalized Gibbs ensemble introduced 
by Rigol et al. in Ref. [H]. 

The Luttinger model (LM) describes a system of inter- 
acting Fermions in one dimension (ID). It was introduced 
by Luttinger j2fli| in 1963, but the correct exact solution 
was found in 1965 by Mattis and Lieb Asymptotic 
forms of one and two-particle correlations in equilibrium 
were obtained by Luther and Peschel 22]. Later, Hal- 
dane |2l [3, [H proposed that this model describes the 
low-energy properties of a fairly broad class of systerns 
in ID known as Tomonaga-Luttinger liquids 23, .26.. 2j|. 

The Hamiltonian of the LM, Hlm = Hq + H2 + H4, 
where Hq = a^'^^FP ■ i''a{p)'4'a{p) ■ is the free-fermion 
Hamiltonian, and the interactions are described by 

9 

H4 = ^^54(9) ■■ MqUai-q) ■■ , (2) 

q.a 

where the Fermi operators {ipaip),ipjj{p')} = SpySa,i3 
(a,/3 — L. R) and anti-commute otherwise. To avoid 
a degenerate ground state, anti-periodic boundary con- 
ditions are chosen: ^l^a{x + L) = —tpaix), {ipa{x) = 
gjsopa: ^p-j ^ ^Yie Fermi field operator and 
sn = —sl = +1 and L the length of the system) so 
that p — 27r(n — \)/L, and n is an integer. The "cur- 
rent" operators Ja[q) — X]q • V'Kp + 9)V'a(p) where 
q = 27rm/L, m being an integer; : . . . ; stands for the 
normal order prescription according to which all creation 
operators are to be found to the left of the annihilation 
operators and expectation values over the ground state 
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of Hq are subtracted. Thus, the above model describes 
a system of fermions interacting via the four Fermion 
terms H2 and H4. Fermions come in two chiralities, R 
standing for right moving and L for left moving parti- 
cles, respectively. The dispersion is linear and therefore 
it is not bounded from below. To define a stable ground 
state for Hq all single-particle levels with p < are filled 
up for both chiralities, which yields a Dirac sea {i.e. an 
"infinite story" hotel) which will be denoted as |0). The 
coupling functions g2{Q) and gi^q) are assumed to be fi- 
nite for q = 0. Moreover, to ensure that the Hilbert 
space of Hlm and Hq remain the same and, in particu- 
lar, that their ground states have a finite overlap at finite 
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as \q\ 
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L, g2{q)/{vF +34(9)) faster than \q 
and 152(g) I <VF+ gi{q) for aU g [H. 

The currents obey a Kac-Moody algebra [3i 
113: [Jo.{q),Jp{q')] = §|<5,+9',o ^a,/3- This fact al- 
lows to introduce, for g 7^ 0, the following operators 
60(g) = -l{2TT/\q\Lf'^[e{q)jR{~q)-e{-q)JL{q)] and 
hl{q) ^l{2^/\q\Lf'^ [e{q)jR{q)-e{-q)JL{~q)l which 
obey the standard algebra of boson ("phonon") oper- 
ators. Moreover, there are two conserved operators 
5N = Nji -\- Nl, i.e. the number operator referred to 
the ground state |0), and the total current J = Nr — Nl, 
where Na = Ja(0). For fermions, the physical states 
obey the selection rule (—1)*^ = (—l)'' . In terms of 
the boson operators 6j(<z), 60(9) the Hamiltonian Hlm is 
quadratic but not diagonal. It can be diagonalized b" 
means of a Bogoliubov ('squeezing') transformation 



i by 



b{q) = coship{q)boiq) +smh f{q)bl{-q), (3) 
b\q) = sinh ip{q) bo{-q) + cosh ip{q) bl{q). (4) 

To render Hlm diagonal, we must choose tSLiih2ip(q) ^ 
g2iq)/[vF + .94(9)]- Thus the Hamiltonian becomes |2J| 
Hlm = Eg^oMlMbHlMl) + ^^^n 5N^/L + 
h-KVj J'^/L, where w(g) = [{vp + gA{q)f - gl{q)V''^, 



s^^', and m 



vj 



being v ~ v{0) and 



vn = ve 
ip = ipiO). 

Let us now consider an interaction quench in the LM. 
Here I consider only the case where the coupling func- 
tions 52(9) and 34(9) are suddenly switched on at i = 0. 
Thus, the initial state of the system will be described by 
a thermal distribution determined by the non-interacting 
Hamiltonian Hq, p{t = 0) = po ~ /Zq, where 

Zq = Tr e~^°/'^ . However, for t > 0, the evolution is 
dictated by the full Hamiltonian Hlm- A more gen- 
eral type of quench corresponds to a sudden switch be- 
tween two different forms of g2(q) and g^iq). Whereas 
the results described below can be generalized to such a 
case, I believe a quench from the non-interacting limit 
is most interesting because the spectrum of Hq con- 
tains free fermions whereas the spectrum of Hlm does 

not m mill 113. Thus, a sudden switch-on of the in- 
teractions describes a time-dependent destruction of the 



characteristic discontinuity of the momentum distribu- 
tion at the Fermi points p — 0. 

Equal time correlations of a given operator 0{x), 

Co{x,t) = (e'^"^*/'*Ot(^)o(o)e-'«^M*/'i)o (5) 
= Trpoe*■f^"^*/'*Ot(x)0(0)e-*■^"^*/^ (6) 

Note that since [Ho,Hlm] 7^ 0, Co{x,t) is exphcitly 
time-dependent. Indeed, in the LM model time depen- 
dence stems from H2, since [Hq,H4] = 0. H2 describes 
scattering between fermions moving in opposite direc- 
tions, and, as shown below, it produces entanglement 
between the excitation modes with q > and q < 0. 
The exact evolution of bo{q) has a fairly simple form: 

bo{q, t) - /(q, i)6o(<z) + .9*(g, t)bl{~q), (7) 
where bo{q,t) = e^Hz^Mt/hfj^^^-^^-zHr^Mt/h^ j^^^^) ^ 
cosv{q)\q\t — i sin v{q)\q\t cosh 2ip{q), and g{q,t) = 
ismv{q)\q\tsmh2Lp{q). Note that this form obeys the 
correct boundary condition bo{q,0) = 60(9)- Entangle- 
ment between modes of opposite q vanishes for ip{q) = 
{i.e. 52(9) = 0) in agreement with the above discussion. 

The evolution of one-particle correlations {i.e. 0{x) = 
'ijja{x)) can be obtained from Eq. |2} and the bosonization 
formula m in 111113: 



ll^a{x) 



(27ra)i/2 



(8) 



being 77^ ^ rjL two different Pauli matrices that 
ensure the anti-commutation of the left and right- 
moving Fermi fields; 4'a{x) ~ Sa^Oa + 2'KxNa/L + 



<i>i{x) + ^a{x), where [N^^ip^p] 



Y^qyo (27r/gL)^/^e^9'^/2e*«^ 6^(9), and a 0+. Setting 
smh2ip{q) = e"''^"'/^ sinh (/j, where i?o is of the 

order of the range of the interactions, and replacing v{q) 
by its q = value [13, simplifies the calculations with- 
out altering the asymptotic form of the correlations. At 
T = 0, for a system of size L, one-body correlations are 
given by the following expression {3lj: 



C^^{x,t>0\L)^G'^°^{x\L) 



Rq 



d{x\L) 



d{x - 2vt\L)d{x + 2vt\L) 
[d{2vt\L)]'^ 



7V2 



(9) 



where d{x\L) = L| sin(7rx/L)|/7r is the cord function, 
G^^\x\L) — i/[2L sin7r(a; + ia)/L] the non-interacing 
correlation function, and 7 = sinh 2ip. The above 
expression is accurate asymptotically, i.e. for d{x ± 
2vt\L),d{x\L),d{2vt\L) :s> Rq. It can be seen that the 
one-particle correlations are periodic in time: C^i^{x,t + 
Tq\L) = C^^{x,t\L) with To = L/2v. This implies 
that the finite-size LM does not relax, which is a con- 
sequence of the (approximately) linear dispersion of the 
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eigenmodes near q = along with the absence of any 
damping mechanisms in the LM (see discussion at the 
end). However, in the thermodynamic limit, L ~* oo, 
and d{x\L) \x\. Therefore, Eq. Q becomes: 



C^^{x,t > 0) 



27r{x + ia) 



{2vty 



7V2 



It is interesting to analyze the above expression in the 
limit where 2vt ^ where it becomes 



iZ{t) 



2-k{x - 



(11) 



being Z(t) — {RQ/2vt)'^ a time-dependent renormal- 
ization constant of the Fermi quasi-particles. Thus for 
short-times the system behaves as a Femi liquid, with a 
singularity at the Fermi points given by Zit), which de- 
creases with time. On the other hand, for 2vt ^ |x| the 
correlation takes a non-Fermi liquid form: 



C^H(it!o < |a;| < 2vt) 



2tt{x 



la) 



Rq 



(12) 



In particular, in the limit t +00 one-particle corre- 
lations relax to the power-law in the right-hand side of 
Eq. ifT^ . Notice that, although C^^(x,i 00) exhibits 
a power-law behavior, the latter is governed by an expo- 
nent that is different from the one that governs asymp- 
totic ground-state correlations [i^, 0| > 7o = 2 sinh^ < 
7^ = sinh^ 2ip for (/? ^ 0. The origin of this new exponent 
will be discussed below. 

The different behavior of C^^ (x, t) for short and long 
times can be understood in terms of a 'light-cone' ef- 
fect ^3: The initial state |0) has higher energy than 
the ground state of i?LM (see discussion further be- 
low). Therefore, it contains long wave-length phonons 
that propagate from time = to time = t along light- 
cones where the role of speed of light is played by v. 
These excitations determine which points retain the same 
type of correlations found in |0) and which points acquire 
new correlations. The latter phenomenon and the over- 
all structure of I|1UJ| bears some resemblance to results 
reported in Ref. ll|. Nevertheless, I have so far failed 
to extend the methods of to the quench in the LM. 
There are two main differences: the initial state in the 
present case is non-critical, and therefore it does not have 
any characteristic (gap) energy scale as the initial states 
considered in [1^ . Secondly, and more importantly, the 
critical exponent found above is different from the bulk 
or boundary exponents of the field operator ipji{x). In- 
deed, this may be an indication that the quench in the 
LM belongs to a different universality class. 

One may think that the relaxation behavior exhib- 
ited by C^j^{x,t) in the thermodynamic limit is be- 
cause the field operator, 'ipii{x), is a non- linear func- 



tion of bo{q) and &J(g). However, the (density) opera- 
tor Jr{x) = dx(j)R{x) /2tt also exhibits relaxation. Set- 
ting 0{x) = Jr{x) in lO, the following is obtained using 
Eq. CH), 



7' 



7 



r}(13) 



2[d{x-2vt)f 2[d{x + 2vt)Y. 

For finite L the density correlation function is again pe- 
riodic in time. However, for i — > cx3, it shows relax- 
ation: Cj^{x,t oo\L) ^ -(1 + 7^)/(47r2x2). This 
form again deviates from the ground state behavior, 
where the prefactor of — l/(47r^x^) is cosh 2iy9 — sinh 2(p — 

It is interesting to find that the above results in the 
t —^ 00 limit can be analytically obtained from the gen- 
eralized Gibbs distribution introduced in Ref. ^2 j which 
is described by the following density matrix: 

P.G = ^e5:^(^)^(«), (14) 



where ZgG = Tre^ Mi)i{q) and [H,I{q)] = [I{q)J{q')] = 
0, that is, a set of independent integrals of motion. Since 
[HLM,'n{q)] = 0, where n{q) — b'^{q)b{q), the phonon 
occupancy operators seem as the most natural choice for 
I{q). The Lagrange multipliers A(g) are obtained from 
the condition |l2l |: 

{niq))t=o = (0|n(q)|0) = {niq))gG = Tr [pgGn{q)] , (15) 

where T — was assumed. Using (|3I4|I . (0|n(g)|0) — 
sinh^ (p{q), which is a non-thermal distribution. However, 
X{q) do not need to obtained explicitly, as it suffices to 
realize that pgc has the same form as the distribution in 
the canonical ensemble with H/T = — J2q '^(9)^1(9)- One 
can also regard pgc as a canonical distribution with a q- 
dependent temperature, T{q) — —hv{q)\q\/ X{q). Using 
this fact along with equations (|3I4|) and JSJ, I find that 

Cfjx) = TrpgG V'li(x)V«(0) = ^]im^C^A^,t), (16) 

C5^(x) = TipgG Jr{x)Jr{0) = ^\im^Cj,{x,t). (17) 

Thus, at least for these simple correlation functions, 
it seems that the generalized Gibbs distribution de- 
scribes the stationary state of the LM after an inter- 
action quench. The reason why the critical exponent 
7^ turns out to be different from the known equilib- 
rium exponents can be thus explained in two different 
ways: mathematically, it is seen that in order to ob- 
tain the evolution of the operator bo{q), Eq one 
has to do and undo the Bogoliubov transformation (|3I4|I . 
However, these transformations do not cancel each other 
exactly (except at t = 0) because of the phase fac- 
tors e^™*^'''l'^l* introduced by the time evolution oper- 
ator. In contrast, in the equilibrium problem, since 
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the expectation value is taken over the ground state 
of Hlm, the Bogohubov transformation is performed 
only once. Physically, in view of the results (|16I17|I . 
the difference in exponent can be regarded a conse- 
quence of the non-equilibrium distribution of phonons 
(0|n(g)|0) = sinh^ ip{q), which is a constant of motion. 
Note as well that {0\nlq)n{q')\0)-{0\n{q)\0){0\n{q')\0) = 
[sinh'* Lp(q) + 2 cosh^ Lp{q) sinh^ fiq)] Sq.q' is non-zero for 
q — q', since |0) is not an eigenstate of Hlm for 52(9) 7^ 0. 

Let us finally consider how the above predictions could 
be experimentally observed. To date, there are no ex- 
act realizations of the LM in nature. However, one can 
exploit the fact that the LM describes the low-energy 
properties of Tomonaga-Luttinger liquids IHIilli^, of 
which several physical realizations in cold atomic gases 
are available |l9l Esl Eof . Let us therefore consider a 
single-species cold Fermi gas confined to ID in a strongly 
anisotropic trap l^^. In a single-species cold Fermi gas, 
the p-wave interaction is naturally negligible. One pos- 
sibility to realize a sudden change of the interaction is 
to use a p-wave Feshbach resonance , which enhances 
the strength of this interaction. Alternatively, one can 
use a ID dipolar Fermi gas, where interactions at long 
distances are described by the potential: 



Vdip(a;) 



1 L>2(l-3cos( 



47reo 



(18) 



D being the dipolar momentum of the atoms and 6 is 
the angle subtended by the direction of the atomic mo- 
tion and an electric field (or magnetic, for magnetic 
dipoles) that polarizes the gas. In the above expres- 
sion Rq is of the order of the transverse size of the 
cloud. The Fourier transform of ((TH|l . 52(9) = 54(9) oc 
Vdipiq) = m\QRo\Kii\qRo\), where A(0) = D^{1 - 
3cos6')/27reoi?o and Ki{x) is the first order modified 
Bessel function. An sudden switch-on of Vdip(g) can be 
realized by deviating the electric field that polarizes the 
gas from the "magic" angle 9m = cos~^(|), for which the 
interaction vanishes (i.e. X{9m) — 0). 

However, the full Hamiltonian for a TLL contains an 
infinite series of terms that spoil the integrability of the 
LM [i^ 0|- Roughly speaking, these stem from the 
non-linearity of the fermion dispersion [32l | and the fact 
that interactions couple right and left-moving modes in 
a way that is highly non-linear in terms of the boson 
fields 4>a{x) (umklapp scattering) In a TLL all 

these deviations are irrelevant in the renormalization- 
group sense, which means that their effect on low-energy 
states is small. Nevertheless, after a sudden change of 
the interaction in the systems described above, high- 
energy excitations will be created that are not described 
by the LM. Exciting many fermions to levels very far 
from the Fermi level where the LM description is not ac- 
curate can be avoided by turning on the interaction to 
a value much smaller than the Fermi energy. On the 
other hand, low energy excitations will survive for longer 



times and, since they dominate the long-time dynam- 
ics, the behavior of the correlations will be described by 
the above results. Thus, if the quench was conducted 
at zero temperature, since the atomic systems are finite, 
an approximately periodic behavior of correlations can 
be expected. However, Fermi gases are usually hard to 
cool down, and a situation where temperature is larger 
than level spacing (i.e. T ^ tt/L) is perfectly reallis- 
tic. In this situation, one should consider correlations 
at finite T, neglecting finite size effects. The latter can 
be obtained from Eq. ^ upon replacing L sin^Trx / L) / it 
by (hvp/TTT) smh^irTx/hvp), etc. Thus relaxation takes 
place because temperature induces a finite correlation 
length in the initial state and therefore correlations de- 
cay exponentially. One-body correlations that can be 
accessed through the momentum distribution, which can 
be measured in a time of light experiment. Thus the 
steady state momentum distribution following a sudden 
switch-on of interactions should differ from the equilib- 
rium distribution at the same tempereature. A more de- 
tailed analysis will be given elsewhere 30] . 
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